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Abstract 

We study theoretically the effect of spin relaxation on the spin transport in a ferromag- 
net/superconductor(FM/SC) tunnel junction. When spin-polarized electrons are injected into the 
SC from the FM, nonequilibrium spin accumulation as well as spin current are created in the range 
of the spin diffusion length in the SC. We find that the spin diffusion length in the superconducting 
state is the same as that in the normal state. We examine a FM/SC/SC double tunnel junction, 
and show that the spin current is detected by the Joule heat generated at the Josephson junction. 
This provides a method to obtain the spin diffusion length by probing the spin current in SC's. 
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Spin transport through magnetic nanostructures has attracted much interest. Using the 
method of tunnehng spectroscopy, Tedrow and Meservey [Q] demonstrated that the current 
through the junction between a ferromagnet (FM) and a superconductor (SC) is spin po- 
larized. Johnson and Silsbee |^ and Jedema et al. Q have observed nonequihbrium 
spin accumulation in a nonmagnetic metal sandwiched by FM's. Suppression of the su- 
perconducting gap due to spin accumulation has been shown experimentally [§],§], and 
theoretically [^. 

SC's are powerful probe for the spin polarization of the current injected from FM's as 
shown in FM/SC tunnel junctions Jll and FM/SC point contacts f^. SC's are also useful for 
exploring how the injected spin-polarized quasiparticles (QP's) are transported, particularly 
the effect of spin relaxation on the spin transport, because the spin-relaxation time and 
the spin diffusion length can be measured precisely in the superconducting state where 
thermal noise effects are extremely small. In addition, the unambiguous description of the 
spin-relaxation effect is possible due to the fact that the spin relaxation is dominated by 
spin-orbit impurity scattering in SC's. 

So far, there have been a number of studies on the spin relaxation time and the spin 
diffusion length in SC's. However, the results are controversial: In a spin coupled resistance 
in permalloy/Nb /permalloy trilayers |Q, it was shown that the spin diffusion length of Nb 
decreases with decreasing temperature in the superconducting state. In contrast, the spin 
relaxation time in SC's was measured by the method of electron spin resonance (ESR) and 
was found to increase with decreasing temperature in SC's [llO| , |Tl| . It was also shown that 
the spin diffusion length in SC's increases with decreasing temperature [|12| by assuming 
that the length is proportional to the square root of the spin relaxation time. Since the 
spin diffusion length and the spin relaxation time are key quantities for the spin transport 
in SC's, it is highly desired to construct a theory of the spin transport and relaxation and 
to solve the controversial issue mentioned above. 

In this paper, we study the spin transport through a FM/SC tunnel junction. The 
spin accumulation and spin current in SC are calculated based on the Boltzmann transport 
theory. It is shown that the spin diffusion length in the superconducting state is equal to 
that in the normal state. We examine a FM/SC1/SC2 double tunnel junction, and show 
that the spin current is detected by the Joule heat generated at the Josephson junction []T3 



which provides information about the spin diffusion length and the spin relaxation time in 
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SC's. 

We first consider a FM/SC tunnel junction. The bias voltage Vt is applied to the tunnel 
junction of resistance i?T- The tunnel barrier is at x = and the current fiows in the x 
direction. The tunnel current is calculated by using the phenomenological tunnel Hamilto- 
nian which describes the transfer of electrons from one electrode to the other. If the SC is 
in the superconducting state, we rewrite the electron operators Oko- in the SC in terms of 
the QP operators 7kcr by using the Bogoliubov transformations, = Mk7kT + 'vl.S'^'^_]^^ and 



a 



^ — -fk5'^7kt +'"k7-ki ; where \u\^\^ = 1 — |fkP = | (1 + (.k/Ey^), S is the operator which 



-ki ~ ^k^ /kT "k /-ki' vviicic |u,k| — -L It^kl — 2 



annihilates a Cooper pair fl^, and i?k = [^k + A^]^/^ is the QP dispersion with ^k and A 
being the one-electron energy relative to the chemical potential of the condensate and the 
superconducting gap, respectively. 

From Fermi's golden rule, the spin-dependent tunnel currents across the FM/SC junction 
are given by 

/TT(^T) = (GTT/e)[Ar-5(0)], (la) 
lTi{V^) = {GTi/e)[Af + S{0)], (lb) 

where Gto- is the tunnel conductance for electrons with spin a when the SC is in the normal 
state, and e the electronic charge. The quantity A/" is the ordinary tunneling term driven 
by^x: f^{VT) = J^^Vs{E)[fo{E-eVT)-foiE)]dE, where VsiE) =Re[\E\/VE^ - A^] 
is the normalized BCS density of states and fo{E) the Fermi distribution function. The 
quantity S{x) is the normalized spin density at position x in the SC; 

5 = l/(2PN)E[/kT-/ki], (2) 

k 

where /ko- = (7ko-7kcr) is the distribution function for a QP with energy E'k and spin a, and 
Pn the density of states in the normal state. In Eq. (1), we neglected the contribution of 
charge imbalance by assuming that the charge diffusion length Aq is much smaller than the 
spin diffusion length A^. It has been reported that A^ ~ 450 /im (Ref. 2) and Ag ~ 10 fim 
(Ref. 15) for Al. The charge current /charge = -^tt + /Ti and the spin current l]^^^ = Jti —Iti 
through the junction are given by 

lL,,, = [M-PSiO)]/ieRT), (3a) 
4pin=[W-5(0)]/(ei?T), (3b) 



where 1/-Rt = Cti + Gtj, and P = (Gtt ~ Gti)/{Gti; + Gti) is the tunnehng spin polar- 
ization. 

Let us examine the effect of spin relaxation on the spin accumulation and the spin current 
in SC. In a steady state, the Boltzmann equation is written as 

vu ■ Vr/k. + k ■ Vk/k. = (5/k./5t),,,tt ' (4) 

where = ^"^Vk-^k = ('Ck/-Ek)fF is the group velocity of QP's and vp the Fermi velocity. 
In the superconducting state, there is no electric field inside SC and thus k = 0. The 
scattering term on the right side of Eq. (^) arose from scattering of QP's by nonmagnetic 
impurities, and is decomposed into the terms due to elastic scattering and spin-flip scattering 

m 

idhJdtlc.^ = -f-^^^ - (5) 

'imp ■^'sf 

where /oo- is the distribution function defined by the average of /ko- with respect to the 
direction of k, Timp = (-Ek/|^k|) Timp ESl, and Tgf = (-Ek/|^k|) ''"if^'' are the elastic and the 



spin-fiip scattering times in the superconducting state, respectively, and t^^^ and r^f' are 
those in the normal state. 

In the FM/SC junction, the physical quantities vary in the x direction and are uniform 
in the yz plane, where Vy/ko- = V^/ko- = 0. From Eqs. (H) and (|), we have 

/ka ~ /o<T - Timp<Vx./oa- (6) 

The spin- dependent current density fiowing in the x direction is calculated as 

poo 

■ V </k. = -2eV^D^^^ / VJo.dE, (7) 



e 

k 



where D^'^^ = ^Vprj^^^ is the diffusion constant in the normal state. 

The spin accumulation at position x in the SC is created by shifting the chemical potential 
of up-spin QP's by 6fi{x) and that of down-spin ones by —6^{x), which is described by taking 
foa{E,x) = fo{E — a6fi{x)). When 6ii is much smaller than A, fo„ is expanded as 

foAE, x) ~ fo{E) - a [dfo{E)/dE] Sfi{x). (8) 

From Eqs. (|^ and (||), the charge current density vanishes: icharge = H + H ~ 0' while the 
spin current density ispin(a;) = — ii is driven by the gradient of 6fi{x), 

^spin(a;) = -4eV^D''-^fo{A)VJ^i{x), (9) 



The divergence of «spin(a;) is given from Eqs. by 
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/o(A)5^(x). 



(10) 



Thus, the chemical potential shift satisfies the equation 



(11) 



where is the spin diffusion length in the SC 




(12) 



In the FM/SC junction, Eq. (|TT]) has a solution of the form Sfi{x) = (5/i(0) exp{—x/Xs), and 
therefore both the spin accumulation and spin current decay exponentially on the length 
scale of Xs- Note that the spin dijfusion length in the superconducting state is the same as that 
in the normal state. This result is understood as follows: In the superconducting state, the 
diffusion constant is D = l^'^'Hrnp = (|'Ck|/-Ek) D^^^ and the spin-flip time Tgf = (-Ek/|^k|) 'Ts{\ 
so that the density of state factor -Ek/|^k| in Xg = \/Dts{ is canceled out, resulting in Eq. (|12|). 

The spin injection experiment has been done to extract the spin diffusion length in 
Nb by using bipolar spin transistors [^]. From the measurement of an excess voltage Vs 
(oc due to spin accumulation, a strong dependence of Vg on temperature (T) was found 
below the superconducting critical temperature Tc. From an analysis of K using the relation 
Vs = Vsoexp{~x/ Xs), As oc (1 — T/Tc)^"^ {n ~ 1/2) was deduced [^f. However, since A^ is 
independent of T as given in Eq. (|l^), A^ in Ref. 3 is not the spin diffusion length, but 
rather the penetration length of the QP evanescent wave into SC due to Andreev reflection 
(AR) [|l^. This is because AR is dominant when SC is in metallic contact with FM's as in 
the experiment of Ref. 3. To measure A^, it is desirable to insert a thin insulating barrier 
between FM and SC for making the QP spin injection predominant. 

Another important quantity for the spin transport is the spin relaxation time Ts of S in 
the superconducting state, which is measured by the ESR experiment. If is introduced 
by the relaxation time approximation {dS /dt)sca.tt = —<S/ts, we find 




(13) 
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For ^ A, Eq. (0) reduces to the result of earlier theories [0. Figure |T] shows the 
temperature dependence of r^. In the normal state (A = 0) above Tc, is equal to T^f\ 
In the superconducting state below T^, increases rapidly with decreasing T and behaves 
similar to ~ {'n-A/2kBT)'^^'^T^^^ at low T. It is worthwhile to note that one cannot use 
Tg in place of t^^'' in Eq. when evaluating A^, because is the relaxation time of 
the macroscopic quantity S while t^^^ is the transport relaxation time of an individual QP 
with particular energy, which makes them different in the superconducting state. In the 
normal state, however, is equal to t^^^ and thus can be used for estimating Ag. Note that 
(5/i(0) ~ (rs/rjj'')(5yu''"-'(0) for eVr <^ A, where 5/i*^"-'(0) is the shift of chemical potential in 
the normal state. 

The above discussions are summarized as follows; (1) The strong T dependence of Vg 
(Ref. 3) is not related to A^ but to the decay length of the evanescent wave in Andreev 
reflection. (2) The ESR experiments [|1^, |ll| are consistent with our theory. (3) The the- 
oretical treatment of A^, in Ref. 12 is not correct because they used the incorrect formula 
A, = \fD^s which differs from Eq. i^. 

In order to investigate the spin diffusion length and the spin current in SC's, we consider 
a FM/SC1/SC2 double tunnel junction. The SCI and SC2 are identical SC's, and their 
thicknesses are d and semi-infinite, respectively. The resistance of the FM/SCl tunnel 
junction and that of the SC1/SC2 Josephson junction (JJ) are i?T and and the voltage 
drops across the junctions are Vt and Vj, respectively. The tunnel current through the J J 
is expressed as 

I = 4'harge (^J) + ^Jl (^j) Sm ^ + /j2 (l^j) COS (14) 

where (p is the phase difference of the gap parameters in SCI and SC2. In Eq. (0), the 
first term describes the QP tunneling, and the second and third terms describe the phase 
coherent (Cooper pair) tunneling. The usual Josephson effect is associated with the snup 
term. Using Fermi's golden rule, we have the spin- dependent QP tunnel current 

1 

/q"p(Kj) = ^ y dEVs{E)Vs{E + eVj) [h{E - a6fi,) - fo{E + eV, - adfi^)] , (15) 

where 5fii{x) is the shift of the chemical potential in the ith SC. The QP charge current 
^charge = ^Ipi^j) + ^<ip(^j) and Spin current 1^^^^ = /^p(Vj) - /j,p(Vj) across the JJ are given 
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by 

Charge = ^ / dEV^{E)V^{E + cV,) [/o(i5; - a5^^{) - /o(i5 + eV, - a5^^^)] , (16) 



spin 



2eR 



^— I dEVs{E)VsiE + eV,)\fo{E - Sfii) - foiE + S^i) 

-/o(^ + eVj - 6fi2) + fo{E + eVj + Sfi^)] ■ (17) 



The phase coherent tunneling terms are obtained as 

A2 ^^ei\E\ - A) ^(A -\E + eVj\) 



'Jl 



'J2 



X 5^ [1 - U\E\ - 5/x,) - /o(|E| + 5/i,)] , (18) 
i=i,2 

^ r rfg Ps(g)Ps(g + e1/.) ^ ^ _ _ _ 

eRj J-oo E{E + eVj) ^ 



where Q{x) is the step function. Equations (|T6D-(|19D are the generahzed formulas for the 



conventional J J [|T^. From Eq. ([TT|) , has the forms b^\[x) m 
SCI and 5iX2{x) = B^e"^^'^^ in SC2, where Bi, B2, and B3 are determined by the boundary 
conditions that the spin currents are continuous at each junction. The results are used to 
calculate the currents through the FM/SC1/SC2 junction. 

In the following we assume that the bias voltage across the J J is zero (Vj = 0). It follows 
from Eqs. ([T6|)-(|T9D that /^harge h2 vanish, whereas Jspj^ oc [5 jjii{d) — 5 iJi2{d)] and Jji ~ Jc, 
Jc being the Josephson critical current. These results indicate that the charge current is 
carried by the Cooper pairs as the dc Josephson current when the bias current is less than 
Jci while the spin current is carried by the QP's as the QP current. 

Figure |] shows the spatial dependence of (— (5/ij) for the up- (down-) spin QP's in 
the SC as well as the pair and QP tunnel currents across the JJ at Vj = 0. The up-spin 
tunnel current across the JJ is driven by the drop [5^i{d) — 5/i2(c?)] in the forward direction, 
while the down-spin one is driven by the same drop in the backward direction. In the SC's, 
the up-spin and down-spin QP's, which are drifted by the slope of the chemical potentials, 
flow in opposite directions to each other, so that the QP's carry only the spin and do not 



carry the charge. This is one of the realizations of spin-charge separation [20 



The most striking feature of the junction is that the spin current across the JJ is accom- 
panied by Joule heating at zero bias voltage (Vj = 0). The power of Joule heating is given 
hy W = [6^i{d) — 6fi2{d)] J/pj^/e, and has the d dependence of the form 

W = P'^a exp{-2d/Xs)/ [1-/3 exp(-2d/A,)]^ (20) 

where 

" e2i?j(l + xi)2(l + 2x2)2' ^ (1 + Xi)(l + 2X2)' 
Here, Xn = (XnfroJVn {n = 1,2), Tq = 2/o(A), and 

r„ = [Vs{E)r (-||) dE, (21) 

with rji = {pj<f\s/ R^A), ri2 = (p^^Xg/RjA), the normal-state resistivity pn of SC, and the 
junction area A [pT| . 

Figure ^ shows the Joule heat as a function of the thickness d of SCI in the case 
where FM is a half metal {P = 1). An efficient generation of W occurs for large values of 
rji, which corresponds to a low area tunnel resistance and/or long A^. It is seen that the 
curves show an exponential decay for c//As > 1; W oc exp{—2d/Xs). At d/Xs = 0.5 and for 
r]i = 0.01, RjA = 10-^ Q cm2, and Aq = 0.39 meV (Al), we obtain W/A = 0.4 mW/cm^ per 
unit area of the JJ, which is large enough to observe experimentally. If W is measured for 
various thickness of SCI at Vj = 0, it provides not only the spin diffusion length As but also 
a direct evidence for the spin current flowing in SC's. Note that our method differs from the 
previous one ; the former probes the spin current and the latter the spin accumulation. 

In conclusion, we have studied the effect of spin relaxation on the spin transport in 
superconductors based on the Boltzmann equation, and shown that the spin diffusion length 
in the superconducting state is equal to that in the normal state. This result resolves the 
controversial issue of the spin diffusion length in the superconducting state. We propose 
a spin-injection device with the Josephson junction to extract information about the spin 
diffusion length and the spin current by measuring Joule heat generated at the Josephson 
junction. 
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FIG. 1: Temperature dependence of the spin relaxation time. The gap Aq is the value of A at 
r = 0. 
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FIG. 2: Spatial variation of the splitting in the chemical potentials of SCI and SC2 in a 
FM/SC1/SC2 tunnel junction. The dashed curves with up and down spins indicate the shifts, 
6fii{x) and —5fii{x), of the up-spin and down-spin quasiparticles (QP's) in the ith SC, respectively, 
and the long dashed line indicates the chemical potential of the Cooper pairs. 



11 




12 



